home *** CD-ROM | disk | FTP | other *** search
/ EuroCD 3 / EuroCD 3.iso / Viewers / aMiPEG_1.0 / src / video.c < prev    next >
C/C++ Source or Header  |  1998-06-24  |  60KB  |  2,233 lines
  1. /*
  2.  *  This file contains C code that implements the video decoder model.
  3.  */
  4.  
  5. #include <cybergraphics/cgxvideo.h>
  6. #include <pragmas/cgxvideo_pragmas.h>
  7. #include <clib/cgxvideo_protos.h>
  8.  
  9. #include <proto/exec.h>
  10. #include <proto/timer.h>
  11. #include <devices/timer.h>
  12.  
  13. #include "video.h"
  14. #include "proto.h"
  15. #include "decoders.h"
  16. #include "util.h"
  17.  
  18. #include <time.h>
  19.  
  20. #define MYRECON
  21.  
  22. /* Declarations of functions. */
  23. //static void init_stat_struct(Statval *);
  24. //static void PrintOneStat(void);
  25. //static void CollectStats(void);
  26. //static unsigned int bitCountRead(void);
  27. //static void StartTime(void);
  28. //static void EndTime(void);
  29. static int ParseSeqHead(VidStream *);
  30. static int ParseGOP(VidStream *);
  31. static int ParsePicture(VidStream *, TimeStamp);
  32. static int ParseSlice(VidStream *);
  33. static int ParseMacroBlock(VidStream *);
  34. static void DoPictureDisplay(VidStream *);
  35.  
  36. #ifndef MYRECON
  37. static void ReconIMBlock(VidStream *, int);
  38. static void ReconPMBlock(VidStream *, int, int, int, int);
  39. static void ReconBMBlock(VidStream *, int, int, int, int);
  40. static void ReconBiMBlock(VidStream *, int, int, int, int, int, int);
  41. #endif
  42.  
  43. static void ProcessSkippedPFrameMBlocks(VidStream *);
  44. static void ProcessSkippedBFrameMBlocks(VidStream *);
  45. static void ReconSkippedBlock16(unsigned char *, unsigned char *, int, int, int, int, int, int, int);
  46. static void ReconSkippedBlock8(unsigned char *, unsigned char *, int, int, int, int, int, int, int);
  47.  
  48. extern int ditherType;
  49.  
  50. /* Macro for returning 1 if num is positive, -1 if negative, 0 if 0. */
  51. //#define Sign(num) ((num > 0) ? 1 : ((num == 0) ? 0 : -1))
  52.  
  53. /* Declare global pointer to vid stream used for current decoding. */
  54.  
  55. VidStream *curVidStream = NULL;
  56.  
  57. /* Array mapping zigzag to array pointer offset. */
  58. int zigzag_direct[64] = {
  59.   0, 1, 8, 16, 9, 2, 3, 10, 17, 24, 32, 25, 18, 11, 4, 5, 12,
  60.   19, 26, 33, 40, 48, 41, 34, 27, 20, 13, 6, 7, 14, 21, 28, 35,
  61.   42, 49, 56, 57, 50, 43, 36, 29, 22, 15, 23, 30, 37, 44, 51,
  62.   58, 59, 52, 45, 38, 31, 39, 46, 53, 60, 61, 54, 47, 55, 62, 63};
  63.  
  64. /* Video rates table (taken from mpeg_play 2.3) */
  65. /* Cheat on Vid rates, round to 30, and use 30 if illegal value 
  66.    Except for 9, where Xing means 15, and given their popularity, we'll
  67.    be nice and do it */
  68. static int VidRateNum[16]={30, 24, 24, 25, 30, 30, 50, 60, 
  69.                          60, 15, 30, 30, 30, 30, 30, 30};
  70.  
  71. /* Initialize P and B skip flags. */
  72.  
  73. static int No_P_Flag = 0;
  74. static int No_B_Flag = 0;
  75.  
  76. /* Max lum, chrom indices for illegal block checking. */
  77.  
  78. static int lmaxx, lmaxy, cmaxx, cmaxy;
  79.  
  80. /* Frame Rate Info */
  81. extern int framerate;
  82.  
  83. /*
  84.   The following accounts for time and size  spent in various parsing acitivites
  85.   if ANALYSIS has been defined.
  86. */
  87. #ifdef ANALYSIS
  88. #include "analysis.h"
  89. #endif
  90.  
  91. double realTimeStart;
  92. int totNumFrames = 0;
  93.  
  94.  
  95. double ReadSysClock(void)
  96. {
  97.     unsigned int clock[2];
  98.  
  99.     (void) timer(clock);
  100.     return clock[0] + clock[1] / 1000000.0;
  101. }
  102.  
  103. void PrintTimeInfo(void)
  104. {
  105.     if (!quietFlag) {
  106.         double spent = ReadSysClock() - realTimeStart;
  107.  
  108.         printf("\nReal Time Spent (After Initializations): %f secs.\n", spent);
  109.         printf("Avg. Frames/Sec: %f\n", ((double) totNumFrames) / spent);
  110.  
  111.         printf("Frame size: %i*%i pixels\n", curVidStream->h_size, curVidStream->v_size);
  112.         if(modeid != 0xffffffff)
  113.             printf("ModeID of display screen: 0x%X\n", modeid);
  114.     }
  115. }
  116.  
  117.  
  118. void ExecuteDisplay(struct vid_stream *vid_stream)
  119. {
  120.   static int rate_deal = -1;
  121.   static int one_frame_time;
  122.  
  123.   static struct timeval tftarget, tfnow;
  124.   int usec, sec;
  125.  
  126.     totNumFrames++;
  127.     if (!quietFlag)
  128.         fprintf (stderr, "%d\r", totNumFrames);
  129.  
  130.  
  131.   /* Do frame rate control */
  132.   switch (rate_deal) {
  133.   case 0:
  134.     break;
  135.   default:
  136.         GetSysTime(&tfnow);
  137.     usec = tftarget.tv_micro - tfnow.tv_micro;
  138.     sec  = tftarget.tv_secs - tfnow.tv_secs;
  139.     if (usec < 0) {
  140.       usec += 1000000;
  141.       sec--;
  142.     }
  143.     
  144.     /* If we're not behind, wait a bit */
  145.     if ((sec >= 0)  &&  usec > 0)
  146.         {
  147.             TimerIO->tr_time.tv_secs = sec;
  148.             TimerIO->tr_time.tv_micro = usec;
  149.  
  150.             DoIO((struct IORequest *)TimerIO);
  151.  
  152.             GetSysTime(&tfnow);
  153.     }
  154.     /* Setup target for next frame */
  155.     tftarget.tv_micro = tfnow.tv_micro + one_frame_time;
  156.     if (tftarget.tv_micro >= 1000000) {
  157.       tftarget.tv_micro -= 1000000;
  158.       tftarget.tv_secs = tfnow.tv_secs + 1;
  159.     } else tftarget.tv_secs = tfnow.tv_secs;
  160.     break;
  161.   case -1:
  162.     switch (framerate) {
  163.     case -1: /* Go with stream Value */
  164.       rate_deal = VidRateNum[vid_stream->picture_rate];
  165.             GetSysTime(&tftarget);
  166.       one_frame_time = 1000000 / rate_deal;
  167.       break;
  168.     case 0: /* as fast as possible */
  169.       rate_deal = 0;
  170.       break;
  171.     default:
  172.       rate_deal = framerate;
  173.             GetSysTime(&tftarget);
  174.       one_frame_time = 1000000 / rate_deal;
  175.       break;
  176.     }
  177.     break;
  178.   }
  179.  
  180.  
  181.     HAM8_draw((char *) vid_stream->current->display, vid_stream->h_size, vid_stream->v_size);
  182. }
  183.  
  184.  
  185. /*
  186.  *--------------------------------------------------------------
  187.  *
  188.  * NewVidStream --
  189.  *
  190.  *    Allocates and initializes a VidStream structure. Takes
  191.  *      as parameter requested size for buffer length.
  192.  *
  193.  * Results:
  194.  *    A pointer to the new VidStream structure.
  195.  *
  196.  * Side effects:
  197.  *      None.
  198.  *
  199.  *--------------------------------------------------------------
  200.  */
  201.  
  202. VidStream * NewVidStream(int bufLength)
  203. {
  204.   int i,j;
  205.   short *intra_quant_matrix, *intra_quant_matrix0, *intra_quant_matrix1;
  206.   VidStream *new;
  207.   static unsigned char default_intra_matrix[64] = {
  208.     8, 16, 19, 22, 26, 27, 29, 34,
  209.     16, 16, 22, 24, 27, 29, 34, 37,
  210.     19, 22, 26, 27, 29, 34, 34, 38,
  211.     22, 22, 26, 27, 29, 34, 37, 40,
  212.     22, 26, 27, 29, 32, 35, 40, 48,
  213.     26, 27, 29, 32, 35, 40, 48, 58,
  214.     26, 27, 29, 34, 38, 46, 56, 69,
  215.     27, 29, 35, 38, 46, 56, 69, 83};
  216.  
  217.   /* Check for legal buffer length. */
  218.  
  219.   if (bufLength < 1024) return NULL;
  220.  
  221.   /* Make buffer length multiple of 4. */
  222.  
  223.   bufLength = (bufLength + 3) >> 2;
  224.  
  225.   /* Allocate memory for new structure. */
  226.  
  227.   new = (VidStream *) malloc(sizeof(VidStream));
  228.  
  229.   /* Initialize pointers to extension and user data. */
  230.  
  231.   new->group.ext_data = new->group.user_data =
  232.     new->picture.extra_info = new->picture.user_data =
  233.     new->picture.ext_data = new->slice.extra_info =
  234.     new->ext_data = new->user_data = NULL;
  235.  
  236.   /* Create default intra matrix and qscale multiplication tables. */
  237.  
  238.   new->intra_quant_matrix_ptr[0] = (short*) malloc(32 * 8*8 * sizeof(short));
  239.   clear64words(intra_quant_matrix = new->intra_quant_matrix_ptr[0]);
  240.   for(j=1; j<32; j++) new->intra_quant_matrix_ptr[j] = (intra_quant_matrix+=64);
  241.  
  242.   intra_quant_matrix0 = intra_quant_matrix = new->intra_quant_matrix_ptr[1];
  243.  
  244.   for(i=0; i<64; i++)
  245.     *intra_quant_matrix++ = default_intra_matrix[zigzag_direct[i]];
  246.  
  247.   for(j=2; j<32; j++)
  248.   {
  249.     intra_quant_matrix1 = new->intra_quant_matrix_ptr[1];
  250.     for(i=0; i<64; i++) *intra_quant_matrix++ = (*intra_quant_matrix0++) + (*intra_quant_matrix1++);
  251.   }
  252.  
  253.   /* Initialize non intra quantization matrix. */
  254.  
  255.   new->non_intra_default = TRUE;
  256.   new->non_intra_quant_matrix_ptr[0] = NULL;
  257.  
  258.   /* Initialize pointers to image spaces. */
  259.  
  260.   new->current = new->past = new->future = NULL;
  261.   for (i = 0; i < RING_BUF_SIZE; i++) new->ring[i] = NULL;
  262.  
  263.   /* Create buffer. */
  264.  
  265.   new->buf_start = (unsigned int *) malloc(bufLength * 4);
  266.  
  267.     if(!new->buf_start)
  268.     {
  269.         fprintf(stderr, "Not enough memory for I/O buffer.\n");
  270.         exit(1);
  271.     }
  272.  
  273.   /*
  274.    * Set max_buf_length to one less than actual length to deal with messy
  275.    * data without proper seq. end codes.
  276.    */
  277.  
  278.   new->max_buf_length = bufLength - 1;
  279.  
  280.   /* Initialize bitstream i/o fields. */
  281.  
  282.   new->bit_offset = 0;
  283.   new->buf_length = 0;
  284.   new->buffer = new->buf_start;
  285.  
  286.   /* display stuff */
  287.  
  288.   new->display_is_initialized = FALSE;
  289.  
  290.   /* Return structure. */
  291.  
  292.   return new;
  293. }
  294.  
  295.  
  296.  
  297. /*
  298.  *--------------------------------------------------------------
  299.  *
  300.  * DestroyVidStream --
  301.  *
  302.  *    Deallocates a VidStream structure.
  303.  *
  304.  * Results:
  305.  *      None.
  306.  *
  307.  * Side effects:
  308.  *    None.
  309.  *
  310.  *--------------------------------------------------------------
  311.  */
  312. void DestroyVidStream(VidStream *astream)
  313. {
  314.   int i;
  315.  
  316.   if (astream->ext_data) free(astream->ext_data);
  317.  
  318.   if (astream->user_data) free(astream->user_data);
  319.  
  320.   if (astream->group.ext_data) free(astream->group.ext_data);
  321.  
  322.   if (astream->group.user_data) free(astream->group.user_data);
  323.  
  324.   if (astream->picture.extra_info) free(astream->picture.extra_info);
  325.  
  326.   if (astream->picture.ext_data) free(astream->picture.ext_data);
  327.  
  328.   if (astream->picture.user_data) free(astream->picture.user_data);
  329.  
  330.   if (astream->slice.extra_info) free(astream->slice.extra_info);
  331.  
  332.   if (astream->buf_start) free(astream->buf_start);
  333.  
  334.   for (i = 0; i < RING_BUF_SIZE; i++) {
  335.     if (astream->ring[i]) {
  336.       DestroyPictImage(astream->ring[i]);
  337.       astream->ring[i] = NULL;
  338.     }
  339.   }
  340.  
  341.   if (!astream->non_intra_default) free(astream->non_intra_quant_matrix_ptr[0]);
  342.  
  343.   if (astream->intra_quant_matrix_ptr[0]) free(astream->intra_quant_matrix_ptr[0]);
  344.  
  345.   free((char *) astream);
  346. }
  347.  
  348.  
  349.  
  350.  
  351. /*
  352.  *--------------------------------------------------------------
  353.  *
  354.  * NewPictImage --
  355.  *
  356.  *    Allocates and initializes a PictImage structure.
  357.  *      The width and height of the image space are passed in
  358.  *      as parameters.
  359.  *
  360.  * Results:
  361.  *    A pointer to the new PictImage structure.
  362.  *
  363.  * Side effects:
  364.  *    None.
  365.  *
  366.  *--------------------------------------------------------------
  367.  */
  368.  
  369. PictImage * NewPictImage(unsigned int width, unsigned int height)
  370. {
  371.   extern int lores, ham6;
  372.   PictImage *new;
  373.   int size = (short)width * (short)height;
  374.  
  375.   /* Allocate memory space for new structure. */
  376.  
  377.   new = (PictImage *) malloc(sizeof(PictImage));
  378.  
  379.   /* Allocate memory for image spaces. */
  380.   if (ditherType == FULL_COLOR_DITHER) {
  381.     new->display = (unsigned char *) malloc(ham6 ? size : (lores ? (size << 1) : (size << 2)));
  382.   } else if (ditherType == CYBERGFX_DITHER || ditherType == CYBERGFXVLAYER_DITHER  || ditherType == CYBERGFXVLAYERGRAY_DITHER) {
  383.         new->display = (unsigned char *) malloc(size * 4);    /* one byte for each of R, G and B  and A :)*/
  384.     } else {
  385.     new->display = (unsigned char *) malloc(size);
  386.     }
  387.  
  388.   new->luminance = (unsigned char *) malloc(size);
  389.   new->Cr = (unsigned char *) malloc(size >> 2);
  390.   new->Cb = (unsigned char *) malloc(size >> 2);
  391.  
  392.   /* Reset locked flag. */
  393.  
  394.   new->locked = 0;
  395.  
  396.   /* Return pointer to new structure. */
  397.  
  398.   return new;
  399. }
  400.  
  401.  
  402.  
  403. /*
  404.  *--------------------------------------------------------------
  405.  *
  406.  * DestroyPictImage --
  407.  *
  408.  *    Deallocates a PictImage structure.
  409.  *
  410.  * Results:
  411.  *      None.
  412.  *
  413.  * Side effects:
  414.  *    None.
  415.  *
  416.  *--------------------------------------------------------------
  417.  */
  418. void DestroyPictImage(PictImage *apictimage)
  419. {
  420.   if (apictimage->luminance) free(apictimage->luminance);
  421.   if (apictimage->Cr) free(apictimage->Cr);
  422.   if (apictimage->Cb) free(apictimage->Cb);
  423.   if (apictimage->display) free(apictimage->display);
  424.   free(apictimage);
  425. }
  426.  
  427. /*
  428.  *--------------------------------------------------------------
  429.  *
  430.  * mpegInitVidRsrc --
  431.  *
  432.  *      reset first to 1
  433.  *
  434.  * gonna check this thingy out !!!!
  435.  *
  436.  *--------------------------------------------------------------
  437.  */
  438. static int first = 1;
  439.  
  440. void mpegInitVidRsrc(void)
  441. {
  442.     first = 1;
  443. }
  444.  
  445.  
  446. /*
  447.  *--------------------------------------------------------------
  448.  *
  449.  * mpegVidRsrc --
  450.  *
  451.  *      Parses bit stream until MB_QUANTUM number of
  452.  *      macroblocks have been decoded or current slice or
  453.  *      picture ends, whichever comes first. If the start
  454.  *      of a frame is encountered, the frame is time stamped
  455.  *      with the value passed in time_stamp. If the value
  456.  *      passed in buffer is not null, the video stream buffer
  457.  *      is set to buffer and the length of the buffer is
  458.  *      expected in value passed in through length. The current
  459.  *      video stream is set to vid_stream. If vid_stream
  460.  *      is passed as NULL, a new VidStream structure is created
  461.  *      and initialized and used as the current video stream.
  462.  *
  463.  * Results:
  464.  *      A Boolean value is returned - 0 is end of picture or error
  465.  *
  466.  * Side effects:
  467.  *      Bit stream is irreversibly parsed. If a picture is completed,
  468.  *      a function is called to display the frame at the correct time.
  469.  *
  470.  *--------------------------------------------------------------
  471.  */
  472.  
  473. int mpegVidRsrc(TimeStamp time_stamp, VidStream *vid_stream)
  474. {
  475.   unsigned int data;
  476.   int i, status;
  477.  
  478.   /* If vid_stream is null, create new VidStream structure.  MR: ????? */
  479.  
  480.     if(!handle_window())
  481.     {
  482.         DestroyVidStream(vid_stream);
  483.         return 0;
  484.     }
  485.  
  486.   if (vid_stream == NULL) return 0;
  487.  
  488.   /*
  489.    * Set global curVid­Stream to vid_stream. Necessary because bit i/o use
  490.    * curVidStream and are not passed vid_stream. Also set global bitstream
  491.    * parameters.
  492.    */
  493.  
  494.   curVidStream = vid_stream;
  495.   bitOffset = curVidStream->bit_offset;
  496.   bufLength = curVidStream->buf_length;
  497.   bitBuffer = curVidStream->buffer;
  498.  
  499.   /*
  500.    * If called for the first time, find start code, make sure it is a
  501.    * sequence start code.
  502.    */
  503.  
  504.   if (first) {
  505.     next_start_code();
  506.     switch(show_bits32(data))
  507.     {
  508.         case SEQ_START_CODE:
  509.         break;
  510.         default:
  511.             fprintf(stderr, "This is not an MPEG stream.\n");
  512.             DestroyVidStream(curVidStream);
  513.             longjmp(env,1);
  514.         break;
  515.     }
  516.     first = 0;
  517.   } 
  518.  
  519.   /*
  520.    * This thing below is a TASTELESS HACK to make system layer streams work.
  521.    * I don't know where the bug lies, but this seems to "fix" it. Somehow the
  522.    * other routines give a load of zeroes sometimes when there should be 0x100.
  523.    * When you unquote the fwrite-stuff in readfile.c to create a non-system
  524.    * layer stream file, it creates an mpeg that doesn't need this...  
  525.    * Still the producing routines are the same! Odd. 
  526.    * -Jussuf the Meddling Master :)
  527.    */
  528.     if(sys_layer)
  529.     {   
  530.         show_bits32(data);
  531.         while(data==0)
  532.         {    
  533.             next_start_code();
  534.             show_bits32(data);
  535.         }
  536.     }
  537.  
  538.   /*
  539.    * Process according to start code (or parse macroblock if not a start code
  540.    * at all.
  541.    */
  542.  
  543.   switch (show_bits32(data)) {
  544.  
  545.   case ISO_11172_END_CODE:
  546.   case SEQ_END_CODE:
  547.  
  548.     /* Display last frame. */
  549.  
  550.     if (vid_stream->future != NULL) {
  551.       vid_stream->current = vid_stream->future;
  552.       ExecuteDisplay(vid_stream);
  553.     }
  554.     
  555.     /* Sequence done. Do the right thing. return False. */
  556.  
  557.     if (!quietFlag) {
  558.       fprintf(stderr, "\nDone!\n");
  559.     }
  560.  
  561. #ifdef ANALYSIS
  562.     PrintAllStats();
  563. #endif
  564.     PrintTimeInfo();
  565.  
  566.     DestroyVidStream(curVidStream);
  567.     return 0;
  568.     break;
  569.  
  570.   case SEQ_START_CODE:
  571.     
  572.     /* Sequence start code. Parse sequence header. */
  573.  
  574.     if (ParseSeqHead(vid_stream) != PARSE_OK) goto error;
  575.  
  576.     /*
  577.      * Return after sequence start code so that application above can use
  578.      * info in header.
  579.      */
  580.  
  581.     goto done;
  582.  
  583.   case GOP_START_CODE:
  584.     /* Group of Pictures start code. Parse gop header. */
  585.  
  586.     if (ParseGOP(vid_stream) != PARSE_OK) goto error;
  587.     goto done;
  588.  
  589.   case PICTURE_START_CODE:
  590.     /* Picture start code. Parse picture header and first slice header. */
  591.  
  592.     status = ParsePicture(vid_stream, time_stamp);
  593.  
  594.     if (status == SKIP_PICTURE) {
  595.       next_start_code();
  596.       fprintf(stderr, "Skipping picture...");
  597.       while ((show_bits32(data)) != PICTURE_START_CODE) {
  598.     if (data == GOP_START_CODE || data == SEQ_END_CODE)
  599.       break;
  600.     flush_bits(24);
  601.     next_start_code();
  602.       }
  603.       fprintf(stderr, "Done.\n");
  604.       goto done;
  605.     } else if (status != PARSE_OK) goto error;
  606.  
  607.  
  608.     if (ParseSlice(vid_stream) != PARSE_OK) goto error;
  609.     break;
  610.  
  611.   case USER_START_CODE:    // there has been one stream with an error in a GOP header
  612.     flush_bits32;
  613.     if (vid_stream->group.user_data) {
  614.       free(vid_stream->group.user_data);
  615.       vid_stream->group.user_data = NULL;
  616.     }
  617.     vid_stream->group.user_data = get_ext_data();
  618.     break;
  619.  
  620.   case EXT_START_CODE:
  621.     flush_bits32;
  622.     if (vid_stream->group.ext_data) {
  623.       free(vid_stream->group.ext_data);
  624.       vid_stream->group.ext_data = NULL;
  625.     }
  626.     vid_stream->group.ext_data = get_ext_data();
  627.     break;
  628.  
  629.   default:
  630.     /* Check for slice start code. */
  631.  
  632.     if ((data >= SLICE_MIN_START_CODE) && (data <= SLICE_MAX_START_CODE))
  633.       if (ParseSlice(vid_stream) != PARSE_OK)    /* Slice start code. Parse slice header. */
  634.     goto error;
  635.  
  636.     break;
  637.   }
  638.  
  639.   /* Parse next MB_QUANTUM macroblocks. */
  640.  
  641.   for (i = 0; i < MB_QUANTUM; i++) {
  642.  
  643.     /* Check to see if actually a startcode and not a macroblock. */
  644.  
  645.    if ((show_bitsn(23,data)) != 0x00000000) {
  646.  
  647.       /* Not start code. Parse Macroblock. */
  648.  
  649.  
  650.   vid_stream->buffer = bitBuffer;
  651.   vid_stream->buf_length = bufLength;
  652.   vid_stream->bit_offset = bitOffset;
  653.  
  654.  
  655.       if (ParseMacroBlock(vid_stream) != PARSE_OK)
  656.     {
  657.        printf("Parsemacroblock error\n");
  658.     goto error;
  659.     }
  660.  
  661. #ifdef ANALYSIS
  662.       if (showmb_flag) {
  663.     DoDitherImage(vid_stream->current->luminance, vid_stream->current->Cr,
  664.               vid_stream->current->Cb, vid_stream->current->display,
  665.               vid_stream->mb_height * 16, vid_stream->mb_width * 16);
  666.     ExecuteDisplay(vid_stream);
  667.       }
  668. #endif
  669.  
  670.     } else {
  671.  
  672.       /* Not macroblock, actually start code. Get start code. */
  673.  
  674.       next_start_code();
  675.       show_bits32(data);
  676.  
  677.       /*
  678.        * If start code is outside range of slice start codes, frame is
  679.        * complete, display frame.
  680.        */
  681.  
  682.       if ((data < SLICE_MIN_START_CODE) || (data > SLICE_MAX_START_CODE)) {
  683.  
  684. #ifdef ANALYSIS
  685.     EndTime();
  686.     stat_a[0].totsize += bitCountRead() - pictureSizeCount;
  687.     if (showEachFlag)
  688.       PrintOneStat();
  689.  
  690.     CollectStats();
  691. #endif
  692.  
  693.         DoPictureDisplay(vid_stream);
  694.       }
  695.       break;
  696.     }
  697.   }
  698.  
  699.   /* Return pointer to video stream structure. */
  700.  
  701.   goto done;
  702.  
  703. error:
  704.   fprintf(stderr, "Error!!!!\n");
  705.   next_start_code();
  706.   goto done;
  707.  
  708. done:
  709.  
  710.   /* Copy global bit i/o variables back into vid_stream. */
  711.  
  712.   vid_stream->buffer = bitBuffer;
  713.   vid_stream->buf_length = bufLength;
  714.   vid_stream->bit_offset = bitOffset;
  715.  
  716.   return 1;
  717. }
  718.  
  719. /*
  720.  *--------------------------------------------------------------
  721.  *
  722.  * ParseSeqHead --
  723.  *
  724.  *      Assumes bit stream is at the begining of the sequence
  725.  *      header start code. Parses off the sequence header.
  726.  *
  727.  * Results:
  728.  *      Fills the vid_stream structure with values derived and
  729.  *      decoded from the sequence header. Allocates the pict image
  730.  *      structures based on the dimensions of the image space
  731.  *      found in the sequence header.
  732.  *
  733.  * Side effects:
  734.  *      Bit stream irreversibly parsed off.
  735.  *
  736.  *--------------------------------------------------------------
  737.  */
  738.  
  739. static int ParseSeqHead(VidStream *vid_stream)
  740. {
  741.   unsigned int data;
  742.   int i;
  743.  
  744.   /* Flush off sequence start code. */
  745.  
  746.   flush_bits32;
  747.  
  748.   /* Get horizontal size of image space. */
  749.  
  750.   get_bits12(vid_stream->h_size);
  751.  
  752.   /* Get vertical size of image space. */
  753.  
  754.   get_bits12(vid_stream->v_size);
  755.  
  756.   /* Calculate macroblock width and height of image space. */
  757.  
  758.   vid_stream->mb_width = (vid_stream->h_size + 15) >> 4;
  759.   vid_stream->mb_height = (vid_stream->v_size + 15) >> 4;
  760.  
  761.   /* Initialize lmaxx, lmaxy, cmaxx, cmaxy. */
  762.  
  763.   lmaxx = (vid_stream->mb_width<<4)-1;
  764.   lmaxy = (vid_stream->mb_height<<4)-1;
  765.   cmaxx = (vid_stream->mb_width<<3)-1;
  766.   cmaxy = (vid_stream->mb_height<<3)-1;
  767.  
  768.   /*
  769.    * Initialize ring buffer of pict images now that dimensions of image space
  770.    * are known.
  771.    */
  772.  
  773.   if (vid_stream->ring[0] == NULL) {
  774.     for (i = 0; i < RING_BUF_SIZE; i++)
  775.       vid_stream->ring[i] = NewPictImage(vid_stream->mb_width << 4, vid_stream->mb_height << 4);
  776.   }
  777.  
  778.   /* Parse of aspect ratio code. */
  779.  
  780.   get_bits4(vid_stream->aspect_ratio);
  781.  
  782.   /* Parse off picture rate code. */
  783.  
  784.   get_bits4(vid_stream->picture_rate);
  785.  
  786.   /* Parse off bit rate. */
  787.  
  788.   get_bits18(vid_stream->bit_rate);
  789.  
  790.   /* Flush marker bit. */
  791.  
  792.   flush_bits(1);
  793.  
  794.   /* Parse off vbv buffer size. */
  795.  
  796.   get_bits10(vid_stream->vbv_buffer_size);
  797.  
  798.   /* Parse off contrained parameter flag. */
  799.  
  800.   get_bits1(vid_stream->const_param_flag);
  801.  
  802.   /*
  803.    * If intra_quant_matrix_flag set, parse off intra quant matrix values.
  804.    */
  805.  
  806.   if (get_bits1(data))
  807.     new_matrix(vid_stream->intra_quant_matrix_ptr[1]);
  808.  
  809.   /*
  810.    * If non intra quant matrix flag set, parse off non intra quant matrix
  811.    * values.
  812.    */
  813.  
  814.   if (get_bits1(data))
  815.   {
  816.     if(vid_stream->non_intra_default)
  817.     {
  818.     unsigned short *intra_quant_matrix;
  819.  
  820.     vid_stream->non_intra_quant_matrix_ptr[0] = (short*) malloc(32 * 8*8 * sizeof(short));
  821.     clear64words(intra_quant_matrix = vid_stream->non_intra_quant_matrix_ptr[0]);
  822.     for(i=1; i<32; i++) vid_stream->non_intra_quant_matrix_ptr[i] = (intra_quant_matrix+=64);
  823.  
  824.     vid_stream->non_intra_default = FALSE;
  825.     }
  826.  
  827.     new_matrix(vid_stream->non_intra_quant_matrix_ptr[1]);
  828.   }
  829.  
  830.   /* Go to next start code. */
  831.  
  832.   next_start_code();
  833.  
  834.   /* If next start code is extension start code, parse off extension data. */
  835.  
  836.   if ((show_bits32(data)) == EXT_START_CODE) {
  837.     flush_bits32;
  838.     if (vid_stream->ext_data) {
  839.       free(vid_stream->ext_data);
  840.       vid_stream->ext_data = NULL;
  841.     }
  842.     vid_stream->ext_data = get_ext_data();
  843.     show_bits32(data);
  844.   }
  845.  
  846.   /* If next start code is user start code, parse off user data. */
  847.  
  848.   if (data == USER_START_CODE) {
  849.     flush_bits32;
  850.     if (vid_stream->user_data) {
  851.       free(vid_stream->user_data);
  852.       vid_stream->user_data = NULL;
  853.     }
  854.     vid_stream->user_data = get_ext_data();
  855.   }
  856.  
  857.   return PARSE_OK;
  858. }
  859.  
  860.  
  861.  
  862. /*
  863.  *--------------------------------------------------------------
  864.  *
  865.  * ParseGOP --
  866.  *
  867.  *      Parses of group of pictures header from bit stream
  868.  *      associated with vid_stream.
  869.  *
  870.  * Results:
  871.  *      Values in gop header placed into video stream structure.
  872.  *
  873.  * Side effects:
  874.  *      Bit stream irreversibly parsed.
  875.  *
  876.  *--------------------------------------------------------------
  877.  */
  878.  
  879. static int ParseGOP(VidStream *vid_stream)
  880. {
  881.   unsigned int data;
  882.  
  883.   /* Flush group of pictures start code. WWWWWWOOOOOOOSSSSSSHHHHH!!! */
  884.  
  885.   flush_bits32;
  886.  
  887.   /* Parse off drop frame flag. */
  888.  
  889.   get_bits1(vid_stream->group.drop_flag);
  890.  
  891.   /* Parse off hour component of time code. */
  892.  
  893.   get_bits5(vid_stream->group.tc_hours);
  894.  
  895.   /* Parse off minute component of time code. */
  896.  
  897.   get_bits6(vid_stream->group.tc_minutes);
  898.  
  899.   /* Flush marker bit. */
  900.  
  901.   flush_bits(1);
  902.  
  903.   /* Parse off second component of time code. */
  904.  
  905.   get_bits6(vid_stream->group.tc_seconds);
  906.  
  907.   /* Parse off picture count component of time code. */
  908.  
  909.   get_bits6(vid_stream->group.tc_pictures);
  910.  
  911.   /* Parse off closed gop and broken link flags. */
  912.  
  913.   get_bits1(vid_stream->group.closed_gop);
  914.   get_bits1(vid_stream->group.broken_link);
  915.  
  916.   /* Goto next start code. */
  917.   next_start_code();
  918.  
  919.   /* If next start code is extension data, parse off extension data. */
  920.  
  921.   if ((show_bits32(data)) == EXT_START_CODE) {
  922.     flush_bits32;
  923.     if (vid_stream->group.ext_data) {
  924.       free(vid_stream->group.ext_data);
  925.       vid_stream->group.ext_data = NULL;
  926.     }
  927.     vid_stream->group.ext_data = get_ext_data();
  928.     show_bits32(data);
  929.   }
  930.  
  931.   /* If next start code is user data, parse off user data. */
  932.  
  933.   if (data == USER_START_CODE) {
  934.     flush_bits32;
  935.     if (vid_stream->group.user_data) {
  936.       free(vid_stream->group.user_data);
  937.       vid_stream->group.user_data = NULL;
  938.     }
  939.     vid_stream->group.user_data = get_ext_data();
  940.   }
  941.  
  942.   return PARSE_OK;
  943. }
  944.  
  945.  
  946.  
  947. /*
  948.  *--------------------------------------------------------------
  949.  *
  950.  * ParsePicture --
  951.  *
  952.  *      Parses picture header. Marks picture to be presented
  953.  *      at particular time given a time stamp.
  954.  *
  955.  * Results:
  956.  *      Values from picture header put into video stream structure.
  957.  *
  958.  * Side effects:
  959.  *      Bit stream irreversibly parsed.
  960.  *
  961.  *--------------------------------------------------------------
  962.  */
  963.  
  964. static int ParsePicture(VidStream *vid_stream, TimeStamp time_stamp)
  965. {
  966.   unsigned int data;
  967.   int i;
  968.  
  969.   /* Flush header start code. */
  970.   flush_bits32;
  971.  
  972.   /* Parse off temporal reference. */
  973.   get_bits10(vid_stream->picture.temp_ref);
  974.  
  975.   /* Parse of picture type. */
  976.   get_bits3(vid_stream->picture.code_type);
  977.  
  978.   if ((vid_stream->picture.code_type == B_TYPE) &&
  979.       (No_B_Flag ||
  980.        (vid_stream->past == NULL) ||
  981.        (vid_stream->future == NULL)))
  982.     return SKIP_PICTURE;
  983.  
  984.   if ((vid_stream->picture.code_type == P_TYPE) &&
  985.       (No_P_Flag || (vid_stream->future == NULL)))
  986.     return SKIP_PICTURE;
  987.  
  988. #ifdef ANALYSIS
  989.   StartTime();
  990.   stat_a[0].frametype = vid_stream->picture.code_type;
  991.   stat_a[0].number = 1;
  992.   stat_a[0].totsize = 45;
  993.   pictureSizeCount = bitCountRead();
  994. #endif
  995.  
  996.   /* Parse off vbv buffer delay value. */
  997.  
  998.   get_bits16(vid_stream->picture.vbv_delay);
  999.  
  1000.   /* If P or B type frame... */
  1001.  
  1002.   if ((vid_stream->picture.code_type == 2) || (vid_stream->picture.code_type == 3)) {
  1003.  
  1004.     /* Parse off forward vector full pixel flag. */
  1005.     get_bits1(vid_stream->picture.full_pel_forw_vector);
  1006.  
  1007.     /* Parse of and decode forw_r_code into forw_r_size and forw_f. */
  1008.  
  1009.     get_bits3(vid_stream->picture.forw_r_size) - 1;
  1010.     vid_stream->picture.forw_f = (1 << vid_stream->picture.forw_r_size);
  1011.   }
  1012.  
  1013.   /* If B type frame... */
  1014.  
  1015.   if (vid_stream->picture.code_type == 3) {
  1016.  
  1017.     /* Parse off back vector full pixel flag. */
  1018.     get_bits1(vid_stream->picture.full_pel_back_vector);
  1019.  
  1020.     /* Parse off and decode back_r_code into back_r_size and back_f. */
  1021.  
  1022.     get_bits3(vid_stream->picture.back_r_size) - 1;
  1023.     vid_stream->picture.back_f = (1 << vid_stream->picture.back_r_size);
  1024.   }
  1025.  
  1026.   /* Get extra bit picture info. */
  1027.  
  1028.   if (vid_stream->picture.extra_info) {
  1029.     free(vid_stream->picture.extra_info);
  1030.     vid_stream->picture.extra_info = NULL;
  1031.   }
  1032.   vid_stream->picture.extra_info = get_extra_bit_info();
  1033.  
  1034.   /* Goto next start code. */
  1035.  
  1036.   next_start_code();
  1037.  
  1038.   /* If start code is extension start code, parse off extension data. */
  1039.  
  1040.   if ((show_bits32(data)) == EXT_START_CODE) {
  1041.     flush_bits32;
  1042.     if (vid_stream->picture.ext_data) {
  1043.       free(vid_stream->picture.ext_data);
  1044.       vid_stream->picture.ext_data = NULL;
  1045.     }
  1046.     vid_stream->picture.ext_data = get_ext_data();
  1047.     show_bits32(data);
  1048.   }
  1049.  
  1050.   /* If start code is user start code, parse off user data. */
  1051.  
  1052.   if (data == USER_START_CODE) {
  1053.     flush_bits32;
  1054.     if (vid_stream->picture.user_data) {
  1055.       free(vid_stream->picture.user_data);
  1056.       vid_stream->picture.user_data = NULL;
  1057.     }
  1058.     vid_stream->picture.user_data = get_ext_data();
  1059.   }
  1060.  
  1061.   /* Find a pict image structure in ring buffer not currently locked. */
  1062.   i = 0;
  1063.   while (vid_stream->ring[i]->locked != 0) {
  1064.     if (++i >= RING_BUF_SIZE) {
  1065.       perror("Fatal error. Ring buffer full.");
  1066.       longjmp(env,1);
  1067.     }
  1068.   }
  1069.  
  1070.   /* Set current pict image structure to the one just found in ring. */
  1071.  
  1072.   vid_stream->current = vid_stream->ring[i];
  1073.  
  1074.   /* Set time stamp. */
  1075.  
  1076.   vid_stream->current->show_time = time_stamp;
  1077.  
  1078.   /* Reset past macroblock address field. */
  1079.  
  1080.   vid_stream->mblock.past_mb_addr = -1;
  1081.  
  1082.   return PARSE_OK;
  1083. }
  1084.  
  1085.  
  1086.  
  1087. /*
  1088.  *--------------------------------------------------------------
  1089.  *
  1090.  * ParseSlice --
  1091.  *
  1092.  *      Parses off slice header.
  1093.  *
  1094.  * Results:
  1095.  *      Values found in slice header put into video stream structure.
  1096.  *
  1097.  * Side effects:
  1098.  *      Bit stream irreversibly parsed.
  1099.  *
  1100.  *--------------------------------------------------------------
  1101.  */
  1102.  
  1103. static int ParseSlice(VidStream *vid_stream)
  1104. {
  1105.   /* Flush slice start code. */
  1106.  
  1107.   flush_bits(24);
  1108.  
  1109.   /* Parse off slice vertical position. */
  1110.  
  1111.   get_bits8(vid_stream->slice.vert_pos);
  1112.  
  1113.   /* Parse off quantization scale. */
  1114.  
  1115.   get_bits5(vid_stream->slice.quant_scale);
  1116.  
  1117.   /* Parse off extra bit slice info. */
  1118.  
  1119.   if (vid_stream->slice.extra_info) {
  1120.     free(vid_stream->slice.extra_info);
  1121.     vid_stream->slice.extra_info = NULL;
  1122.   }
  1123.   vid_stream->slice.extra_info = get_extra_bit_info();
  1124.  
  1125.   /* Reset past intrablock address. */
  1126.  
  1127.   vid_stream->mblock.past_intra_addr = -2;
  1128.  
  1129.   /* Reset previous recon motion vectors. */
  1130.  
  1131.   vid_stream->mblock.recon_right_for_prev = 0;
  1132.   vid_stream->mblock.recon_down_for_prev = 0;
  1133.   vid_stream->mblock.recon_right_back_prev = 0;
  1134.   vid_stream->mblock.recon_down_back_prev = 0;
  1135.  
  1136.   /* Reset macroblock address. */
  1137.  
  1138.   vid_stream->mblock.mb_address = ((short)(vid_stream->slice.vert_pos - 1) * (short)vid_stream->mb_width) - 1;
  1139.  
  1140.   /* Reset past dct dc y, cr, and cb values. */
  1141.  
  1142.   vid_stream->block.dct_dc_y_past = 1024;
  1143.   vid_stream->block.dct_dc_cr_past = 1024;
  1144.   vid_stream->block.dct_dc_cb_past = 1024;
  1145.  
  1146.   return PARSE_OK;
  1147. }
  1148.  
  1149.  
  1150.  
  1151. /*
  1152.  *--------------------------------------------------------------
  1153.  *
  1154.  * ParseMacroBlock --
  1155.  *
  1156.  *      Parseoff macroblock. Reconstructs DCT values. Applies
  1157.  *      inverse DCT, reconstructs motion vectors, calculates and
  1158.  *      set pixel values for macroblock in current pict image
  1159.  *      structure.
  1160.  *
  1161.  * Results:
  1162.  *      Here's where everything really happens. Welcome to the
  1163.  *      heart of darkness.
  1164.  *
  1165.  * Side effects:
  1166.  *      Bit stream irreversibly parsed off.
  1167.  *
  1168.  *--------------------------------------------------------------
  1169.  */
  1170.  
  1171. static int ParseMacroBlock(VidStream *vid_stream)
  1172. {
  1173.     int addr_incr, mask, i;
  1174.     int recon_right_for, recon_down_for, recon_right_back, recon_down_back, zero_block_flag;
  1175.     int mb_data;
  1176.  
  1177. #ifdef ANALYSIS
  1178.     mbSizeCount = bitCountRead();
  1179. #endif
  1180.  
  1181.     /*
  1182.      * Parse off macroblock address increment and add to macroblock address.
  1183.      */
  1184.     do {
  1185.         DecodeMBAddrInc(addr_incr);
  1186.         if (addr_incr == MB_ESCAPE) {
  1187.             vid_stream->mblock.mb_address += 33;
  1188.             addr_incr = MB_STUFFING;
  1189.         }
  1190.     } while (addr_incr == MB_STUFFING);
  1191.     vid_stream->mblock.mb_address += addr_incr;
  1192.  
  1193.     if (vid_stream->mblock.mb_address > ((short)vid_stream->mb_height * (short)vid_stream->mb_width - 1))
  1194.         return SKIP_TO_START_CODE;
  1195.  
  1196.     /*
  1197.      * If macroblocks have been skipped, process skipped macroblocks.
  1198.      */
  1199.  
  1200.     if (vid_stream->mblock.mb_address - vid_stream->mblock.past_mb_addr > 1) {
  1201.         if (vid_stream->picture.code_type == P_TYPE)
  1202.             ProcessSkippedPFrameMBlocks(vid_stream);
  1203.         else if (vid_stream->picture.code_type == B_TYPE)
  1204.             ProcessSkippedBFrameMBlocks(vid_stream);
  1205.     }
  1206.  
  1207.     /* Set past macroblock address to current macroblock address. */
  1208.     vid_stream->mblock.past_mb_addr = vid_stream->mblock.mb_address;
  1209.  
  1210.     /* Based on picture type decode macroblock type. */
  1211.     switch (vid_stream->picture.code_type) {
  1212.         case I_TYPE:
  1213.             DecodeMBTypeI(mb_data, vid_stream->mblock.mb_intra);
  1214.             break;
  1215.         case P_TYPE:
  1216.             DecodeMBTypeP(mb_data, vid_stream->mblock.mb_intra);
  1217.             break;
  1218.         case B_TYPE:
  1219.             DecodeMBTypeB(mb_data, vid_stream->mblock.mb_intra);
  1220.             break;
  1221.     }
  1222.  
  1223.     /* If quantization flag set, parse off new quantization scale. */
  1224.     if (mb_data&MB_quant) {
  1225.         get_bits5(vid_stream->slice.quant_scale);
  1226.     }
  1227.  
  1228.   /* If forward motion vectors exist... */
  1229.   if (mb_data&MB_motion_forw) {
  1230.  
  1231.     /* Parse off and decode horizontal forward motion vector. */
  1232.     DecodeMotionVectors(vid_stream->mblock.motion_h_forw_code);
  1233.  
  1234.     /* If horiz. forward r data exists, parse off. */
  1235.  
  1236.     if ((vid_stream->picture.forw_f != 1) && (vid_stream->mblock.motion_h_forw_code != 0)) {
  1237.       get_bitsn(vid_stream->picture.forw_r_size, vid_stream->mblock.motion_h_forw_r);
  1238.     }
  1239.     /* Parse off and decode vertical forward motion vector. */
  1240.     DecodeMotionVectors(vid_stream->mblock.motion_v_forw_code);
  1241.  
  1242.     /* If vert. forw. r data exists, parse off. */
  1243.  
  1244.     if ((vid_stream->picture.forw_f != 1) && (vid_stream->mblock.motion_v_forw_code != 0)) {
  1245.       get_bitsn(vid_stream->picture.forw_r_size, vid_stream->mblock.motion_v_forw_r);
  1246.     }
  1247.   }
  1248.   /* If back motion vectors exist... */
  1249.   if (mb_data&MB_motion_back) {
  1250.  
  1251.     /* Parse off and decode horiz. back motion vector. */
  1252.     DecodeMotionVectors(vid_stream->mblock.motion_h_back_code);
  1253.  
  1254.     /* If horiz. back r data exists, parse off. */
  1255.  
  1256.     if ((vid_stream->picture.back_f != 1) && (vid_stream->mblock.motion_h_back_code != 0)) {
  1257.       get_bitsn(vid_stream->picture.back_r_size, vid_stream->mblock.motion_h_back_r);
  1258.     }
  1259.     /* Parse off and decode vert. back motion vector. */
  1260.     DecodeMotionVectors(vid_stream->mblock.motion_v_back_code);
  1261.  
  1262.     /* If vert. back r data exists, parse off. */
  1263.  
  1264.     if ((vid_stream->picture.back_f != 1) && (vid_stream->mblock.motion_v_back_code != 0)) {
  1265.       get_bitsn(vid_stream->picture.back_r_size, vid_stream->mblock.motion_v_back_r);
  1266.     }
  1267.   }
  1268. #ifdef ANALYSIS
  1269.     if (vid_stream->mblock.mb_intra) {
  1270.         stat_a[0].i_mbnum++;
  1271.         mbCBPPtr = stat_a[0].i_mbcbp;
  1272.         mbCoeffPtr = stat_a[0].i_mbcoeff;
  1273.         mbSizePtr = &(stat_a[0].i_mbsize);
  1274.     } else if ((mb_data&(MB_motion_back | MB_motion_forw)) == (MB_motion_back | MB_motion_forw)) {
  1275.         stat_a[0].bi_mbnum++;
  1276.         mbCBPPtr = stat_a[0].bi_mbcbp;
  1277.         mbCoeffPtr = stat_a[0].bi_mbcoeff;
  1278.         mbSizePtr = &(stat_a[0].bi_mbsize);
  1279.     } else if (mb_data&MB_motion_back) {
  1280.         stat_a[0].b_mbnum++;
  1281.         mbCBPPtr = stat_a[0].b_mbcbp;
  1282.         mbCoeffPtr = stat_a[0].b_mbcoeff;
  1283.         mbSizePtr = &(stat_a[0].b_mbsize);
  1284.     } else {
  1285.         stat_a[0].p_mbnum++;
  1286.         mbCBPPtr = stat_a[0].p_mbcbp;
  1287.         mbCoeffPtr = stat_a[0].p_mbcoeff;
  1288.         mbSizePtr = &(stat_a[0].p_mbsize);
  1289.     }
  1290. #endif
  1291.  
  1292.     /* If mblock pattern flag set, parse and decode CBP (code block pattern). */
  1293.     if (mb_data&MB_pattern) {
  1294.         DecodeCBP(vid_stream->mblock.cbp);
  1295.     } else /* Otherwise, set CBP to zero. */
  1296.         vid_stream->mblock.cbp = 0;
  1297.  
  1298. #ifdef ANALYSIS
  1299.     mbCBPPtr[vid_stream->mblock.cbp]++;
  1300. #endif
  1301.  
  1302.   /* Reconstruct motion vectors depending on picture type. */
  1303.   if (vid_stream->picture.code_type == P_TYPE) {
  1304.  
  1305.     /*
  1306.      * If no forw motion vectors, reset previous and current vectors to 0.
  1307.      */
  1308.  
  1309.     if (!(mb_data&MB_motion_forw)) {
  1310.       recon_right_for = 0;
  1311.       recon_down_for = 0;
  1312.       vid_stream->mblock.recon_right_for_prev = 0;
  1313.       vid_stream->mblock.recon_down_for_prev = 0;
  1314.     }
  1315.     /*
  1316.      * Otherwise, compute new forw motion vectors. Reset previous vectors to
  1317.      * current vectors.
  1318.      */
  1319.  
  1320.     else {
  1321.       ComputeForwVector(&recon_right_for, &recon_down_for);
  1322.     }
  1323.   }
  1324.   if (vid_stream->picture.code_type == B_TYPE) {
  1325.  
  1326.     /* Reset prev. and current vectors to zero if mblock is intracoded. */
  1327.  
  1328.     if (vid_stream->mblock.mb_intra) {
  1329.       vid_stream->mblock.recon_right_for_prev = 0;
  1330.       vid_stream->mblock.recon_down_for_prev = 0;
  1331.       vid_stream->mblock.recon_right_back_prev = 0;
  1332.       vid_stream->mblock.recon_down_back_prev = 0;
  1333.     } else {
  1334.  
  1335.       /* If no forw vectors, current vectors equal prev. vectors. */
  1336.  
  1337.       if (!(mb_data&MB_motion_forw)) {
  1338.     recon_right_for = vid_stream->mblock.recon_right_for_prev;
  1339.     recon_down_for = vid_stream->mblock.recon_down_for_prev;
  1340.       }
  1341.       /*
  1342.        * Otherwise compute forw. vectors. Reset prev vectors to new values.
  1343.        */
  1344.  
  1345.       else {
  1346.         ComputeForwVector(&recon_right_for, &recon_down_for);
  1347.       }
  1348.  
  1349.       /* If no back vectors, set back vectors to prev back vectors. */
  1350.  
  1351.       if (!(mb_data&MB_motion_back)) {
  1352.     recon_right_back = vid_stream->mblock.recon_right_back_prev;
  1353.     recon_down_back = vid_stream->mblock.recon_down_back_prev;
  1354.       }
  1355.       /* Otherwise compute new vectors and reset prev. back vectors. */
  1356.  
  1357.       else {
  1358.     ComputeBackVector(&recon_right_back, &recon_down_back);
  1359.       }
  1360.  
  1361.       /*
  1362.        * Store vector existance flags in structure for possible skipped
  1363.        * macroblocks to follow.
  1364.        */
  1365.  
  1366.       if(mb_data<0)
  1367.     vid_stream->mblock.bpict_past_forw = vid_stream->mblock.bpict_past_back = -1;
  1368.       else {
  1369.     vid_stream->mblock.bpict_past_forw = mb_data&MB_motion_forw;
  1370.     vid_stream->mblock.bpict_past_back = mb_data&MB_motion_back;
  1371.       }
  1372.  
  1373.     }
  1374.   }
  1375.  
  1376.   /* For each possible block in macroblock. */
  1377.   if (ditherType == GRAY_DITHER || ditherType == CYBERGFXGRAY_DITHER  || ditherType == CYBERGFXVLAYERGRAY_DITHER) {
  1378.     for (mask = 32, i = 0; i < 4; mask >>= 1, i++) {
  1379.  
  1380.       /* If block exists... */
  1381.       if ((vid_stream->mblock.mb_intra) || (vid_stream->mblock.cbp & mask)) {
  1382.     zero_block_flag = 0;
  1383.     ParseReconBlock(i);
  1384.       } else
  1385.     zero_block_flag = 1;
  1386.  
  1387.       /* If macroblock is intra coded... */
  1388.       if (vid_stream->mblock.mb_intra) {
  1389.     ReconIMBlock(vid_stream, i);
  1390.       } else if ((mb_data&(MB_motion_forw | MB_motion_back)) == (MB_motion_forw | MB_motion_back)) {
  1391.     ReconBiMBlock(vid_stream, i, recon_right_for, recon_down_for,
  1392.               recon_right_back, recon_down_back, zero_block_flag);
  1393.       } else if ((mb_data&MB_motion_forw) || (vid_stream->picture.code_type == P_TYPE)) {
  1394.     ReconPMBlock(vid_stream, i, recon_right_for, recon_down_for,
  1395.              zero_block_flag);
  1396.       } else if (mb_data&MB_motion_back) {
  1397.     ReconBMBlock(vid_stream, i, recon_right_back, recon_down_back,
  1398.              zero_block_flag);
  1399.       }
  1400.     }
  1401.     /* Kill the Chrominace blocks... */
  1402.     if ((vid_stream->mblock.mb_intra) || (vid_stream->mblock.cbp & 0x2)) ParseAwayBlock(4);
  1403.     if ((vid_stream->mblock.mb_intra) || (vid_stream->mblock.cbp & 0x1)) ParseAwayBlock(5);
  1404.   } else {
  1405.       for (mask = 32, i = 0; i < 6; mask >>= 1, i++) {
  1406.     
  1407.     /* If block exists... */
  1408.     if ((vid_stream->mblock.mb_intra) || (vid_stream->mblock.cbp & mask)) {
  1409.       zero_block_flag = 0;
  1410.       ParseReconBlock(i);
  1411.     } else
  1412.       zero_block_flag = 1;
  1413.     
  1414.     /* If macroblock is intra coded... */
  1415.     if (vid_stream->mblock.mb_intra) {
  1416.       ReconIMBlock(vid_stream, i);
  1417.     } else if ((mb_data&(MB_motion_forw | MB_motion_back)) == (MB_motion_forw | MB_motion_back)) {
  1418.       ReconBiMBlock(vid_stream, i, recon_right_for, recon_down_for,
  1419.             recon_right_back, recon_down_back, zero_block_flag);
  1420.     } else if ((mb_data&MB_motion_forw) || (vid_stream->picture.code_type == P_TYPE)) {
  1421.       ReconPMBlock(vid_stream, i, recon_right_for, recon_down_for,
  1422.                zero_block_flag);
  1423.     } else if (mb_data&MB_motion_back) {
  1424.       ReconBMBlock(vid_stream, i, recon_right_back, recon_down_back,
  1425.                zero_block_flag);
  1426.     }
  1427.       }
  1428.   }
  1429.  
  1430.   /* If D Type picture, flush marker bit. */
  1431.   if (vid_stream->picture.code_type == 4) flush_bits(1);
  1432.  
  1433.   /* If macroblock was intracoded, set macroblock past intra address. */
  1434.   if (vid_stream->mblock.mb_intra) vid_stream->mblock.past_intra_addr = vid_stream->mblock.mb_address;
  1435.  
  1436. #ifdef ANALYSIS
  1437.   *mbSizePtr += bitCountRead() - mbSizeCount;
  1438. #endif
  1439.  
  1440.   return PARSE_OK;
  1441. }
  1442.  
  1443. #ifndef MYRECON
  1444.  
  1445. /*MR
  1446.  *--------------------------------------------------------------
  1447.  *
  1448.  * ReconIMBlock --
  1449.  *
  1450.  *    Reconstructs intra coded macroblock.
  1451.  *
  1452.  * Results:
  1453.  *    None.
  1454.  *
  1455.  * Side effects:
  1456.  *    None.
  1457.  *
  1458.  *--------------------------------------------------------------
  1459.  */
  1460.  
  1461. static void ReconIMBlock(VidStream *vid_stream, int bnum)
  1462. {
  1463.   int row, col, row_size;
  1464.   unsigned char *dest;
  1465.  
  1466.   row_size = vid_stream->mb_width;
  1467.   /* Calculate macroblock row and column from address. */
  1468.   row = vid_stream->mblock.mb_address / row_size;
  1469.   col = vid_stream->mblock.mb_address % row_size;
  1470.  
  1471.   /* If block is luminance block... */
  1472.  
  1473.   if (bnum < 4) {
  1474.  
  1475.     /* Calculate row and col values for upper left pixel of block. */
  1476.     row <<= 4;
  1477.     col <<= 4;
  1478.     if (bnum > 1) row += 8;
  1479.     if (bnum & 1) col += 8;
  1480.  
  1481.     /* Set dest to luminance plane of current pict image. */
  1482.     dest = vid_stream->current->luminance;
  1483.  
  1484.     /* Establish row size. */
  1485.     row_size *= 16;
  1486.   }
  1487.   /* Otherwise if block is Cr block... */
  1488.  
  1489.   else if (bnum == 4) {
  1490.  
  1491.     /* Set dest to Cr plane of current pict image. */
  1492.     dest = vid_stream->current->Cr;
  1493.  
  1494.     /* Establish row size. */
  1495.     row_size <<= 3;
  1496.  
  1497.     /* Calculate row,col for upper left pixel of block. */
  1498.     row <<= 3;
  1499.     col <<= 3;
  1500.   }
  1501.   /* Otherwise block is Cb block, and ... */
  1502.  
  1503.   else {
  1504.  
  1505.     /* Set dest to Cb plane of current pict image. */
  1506.     dest = vid_stream->current->Cb;
  1507.  
  1508.     /* Establish row size. */
  1509.     row_size <<= 3;
  1510.  
  1511.     /* Calculate row,col for upper left pixel value of block. */
  1512.     row <<= 3;
  1513.     col <<= 3;
  1514.   }
  1515.  
  1516.   /*
  1517.    * For each pixel in block, set to cropped reconstructed value from inverse dct.
  1518.    */
  1519.   IM_reconstruct(dest + row * row_size + col, &vid_stream->block.dct_recon[0][0], row_size);
  1520. }
  1521.  
  1522.  
  1523. /* MR
  1524.  *--------------------------------------------------------------
  1525.  *
  1526.  * ReconPMBlock --
  1527.  *
  1528.  *    Reconstructs forward predicted macroblocks.
  1529.  *
  1530.  * Results:
  1531.  *      None.
  1532.  *
  1533.  * Side effects:
  1534.  *      None.
  1535.  *
  1536.  *--------------------------------------------------------------
  1537.  */
  1538. static void ReconPMBlock(VidStream *vid_stream, int bnum, int recon_right_for, int recon_down_for, int zflag)
  1539. {
  1540.   int row, col, row_size;
  1541.   unsigned char *dest, *past, *rindex2;
  1542.   int right_half_for, down_half_for;
  1543.  
  1544.   row_size = vid_stream->mb_width;
  1545.   /* Calculate macroblock row and column from address. */
  1546.   row = vid_stream->mblock.mb_address / row_size;
  1547.   col = vid_stream->mblock.mb_address % row_size;
  1548.  
  1549.   if (bnum < 4) {
  1550.     /* Calculate right_for, down_for motion vectors. */
  1551.  
  1552.     right_half_for = recon_right_for & 0x1;
  1553.     down_half_for = recon_down_for & 0x1;
  1554.     recon_right_for >>= 1;
  1555.     recon_down_for >>= 1;
  1556.  
  1557.     /* Set dest to luminance plane of current pict image. */
  1558.     dest = vid_stream->current->luminance;
  1559.  
  1560.     if (vid_stream->picture.code_type == B_TYPE) {
  1561.       if (vid_stream->past) past = vid_stream->past->luminance;
  1562.     } else {
  1563.       /* Set predicitive frame to current future frame. */
  1564.       if (vid_stream->future) past = vid_stream->future->luminance;
  1565.     }
  1566.  
  1567.     /* Establish row size. */
  1568.     row_size <<= 4;
  1569.  
  1570.     /* Calculate row,col of upper left pixel in block. */
  1571.     row <<= 4;
  1572.     col <<= 4;
  1573.     if (bnum > 1) row += 8;
  1574.     if (bnum & 1) col += 8;
  1575.  
  1576.   }
  1577.   /* Otherwise, block is NOT luminance block, ... */
  1578.   else {
  1579.     /* Construct motion vectors. */
  1580.     right_half_for = (recon_right_for & 0x2)>>1;
  1581.     down_half_for = recon_down_for & 0x2;
  1582.     recon_right_for >>= 2;
  1583.     recon_down_for >>= 2;
  1584.  
  1585.     /* Establish row size. */
  1586.     row_size <<= 3;
  1587.  
  1588.     /* Calculate row,col of upper left pixel in block. */
  1589.     row <<= 3;
  1590.     col <<= 3;
  1591.  
  1592.     /* If block is Cr block... */
  1593.  
  1594.     if (bnum == 4) {
  1595.       /* Set dest to Cr plane of current pict image. */
  1596.       dest = vid_stream->current->Cr;
  1597.  
  1598.       if (vid_stream->picture.code_type == B_TYPE) {
  1599.     if (vid_stream->past) past = vid_stream->past->Cr;
  1600.       } else {
  1601.     if (vid_stream->future) past = vid_stream->future->Cr;
  1602.       }
  1603.     }
  1604.     /* Otherwise, block is Cb block... */
  1605.     else {
  1606.       /* Set dest to Cb plane of current pict image. */
  1607.       dest = vid_stream->current->Cb;
  1608.  
  1609.       if (vid_stream->picture.code_type == B_TYPE) {
  1610.     if (vid_stream->past) past = vid_stream->past->Cb;
  1611.       } else {
  1612.     if (vid_stream->future) past = vid_stream->future->Cb;
  1613.       }
  1614.     }
  1615.   }
  1616.  
  1617.   /* For each pixel in block... */
  1618.     dest += ((short)row * (short)row_size) + col;
  1619.     past += (short)(row + recon_down_for) * (short)row_size + col + recon_right_for;
  1620.  
  1621.     /*
  1622.      * Calculate predictive pixel value based on motion vectors and copy to
  1623.      * dest plane.
  1624.      */
  1625.  
  1626.     if ((!down_half_for) && (!right_half_for)) {
  1627.       if (!zflag)
  1628.     BMcm_reconstruct(dest, past, &(vid_stream->block.dct_recon[0][0]), row_size);
  1629.       else
  1630.     BM_reconstruct(dest, past, row_size);
  1631.     } else {
  1632.       rindex2 = past + right_half_for + (down_half_for?row_size:0);
  1633.       if (!zflag)
  1634.     BIMcm_reconstruct(dest, past, rindex2, &(vid_stream->block.dct_recon[0][0]), row_size);
  1635.       else
  1636.     BIM_reconstruct(dest, past, rindex2, row_size);
  1637.     }
  1638.  
  1639. }
  1640.  
  1641.  
  1642. /* MR
  1643.  *--------------------------------------------------------------
  1644.  *
  1645.  * ReconBMBlock --
  1646.  *
  1647.  *    Reconstructs back predicted macroblocks.
  1648.  *
  1649.  * Results:
  1650.  *      None.
  1651.  *
  1652.  * Side effects:
  1653.  *      None.
  1654.  *
  1655.  *--------------------------------------------------------------
  1656.  */
  1657.  
  1658. static void ReconBMBlock(VidStream *vid_stream, int bnum, int recon_right_back, int recon_down_back, int zflag)
  1659. {
  1660.   int row, col, row_size;
  1661.   unsigned char *dest, *future, *rindex2;
  1662.   int right_half_back, down_half_back;
  1663.  
  1664.   row_size = vid_stream->mb_width;
  1665.   /* Calculate macroblock row and column from address. */
  1666.   row = vid_stream->mblock.mb_address / row_size;
  1667.   col = vid_stream->mblock.mb_address % row_size;
  1668.  
  1669.   /* If block is luminance block... */
  1670.  
  1671.   if (bnum < 4) {
  1672.     /* Calculate right_back, down_back motion vectors. */
  1673.     right_half_back = recon_right_back & 0x1;
  1674.     down_half_back = recon_down_back & 0x1;
  1675.     recon_right_back >>= 1;
  1676.     recon_down_back >>= 1;
  1677.  
  1678.     /* Set dest to luminance plane of current pict image. */
  1679.     dest = vid_stream->current->luminance;
  1680.  
  1681.     /* If future frame exists, set future to luminance plane of future frame. */
  1682.     if (vid_stream->future) future = vid_stream->future->luminance;
  1683.  
  1684.     /* Establish row size. */
  1685.     row_size <<= 4;
  1686.  
  1687.     /* Calculate row,col of upper left pixel in block. */
  1688.     row <<= 4;
  1689.     col <<= 4;
  1690.     if (bnum > 1) row += 8;
  1691.     if (bnum & 1) col += 8;
  1692.  
  1693.   }
  1694.   /* Otherwise, block is NOT luminance block, ... */
  1695.   else {
  1696.     /* Construct motion vectors. */
  1697.     right_half_back = (recon_right_back & 0x2)>>1;    /* used as integer down there */
  1698.     down_half_back = recon_down_back & 0x2;
  1699.     recon_right_back >>= 2;
  1700.     recon_down_back >>= 2;
  1701.  
  1702.     /* Establish row size. */
  1703.     row_size <<= 3;
  1704.  
  1705.     /* Calculate row,col of upper left pixel in block. */
  1706.     row <<= 3;
  1707.     col <<= 3;
  1708.  
  1709.     /* If block is Cr block... */
  1710.  
  1711.     if (bnum == 4) {
  1712.       /* Set dest to Cr plane of current pict image. */
  1713.       dest = vid_stream->current->Cr;
  1714.  
  1715.       /* If future frame exists, set future to Cr plane of future image. */
  1716.       if (vid_stream->future) future = vid_stream->future->Cr;
  1717.     }
  1718.     /* Otherwise, block is Cb block... */
  1719.     else {
  1720.       /* Set dest to Cb plane of current pict image. */
  1721.       dest = vid_stream->current->Cb;
  1722.  
  1723.       /* If future frame exists, set future to Cb plane of future frame. */
  1724.       if (vid_stream->future) future = vid_stream->future->Cb;
  1725.     }
  1726.   }
  1727.  
  1728.   /* For each pixel in block do... */
  1729.     dest += ((short)row * (short)row_size) + col;
  1730.     future += (short)(row + recon_down_back) * (short)row_size + col + recon_right_back;
  1731.  
  1732.     if ((!right_half_back) && (!down_half_back)) {
  1733.       if (!zflag)
  1734.     BMcm_reconstruct(dest, future, &(vid_stream->block.dct_recon[0][0]), row_size);
  1735.       else
  1736.     BM_reconstruct(dest, future, row_size);
  1737.     } else {
  1738.       rindex2 = future + right_half_back + (down_half_back?row_size:0);
  1739.       if (!zflag)
  1740.     BIMcm_reconstruct(dest, future, rindex2, &(vid_stream->block.dct_recon[0][0]), row_size);
  1741.       else
  1742.     BIM_reconstruct(dest, future, rindex2, row_size);
  1743.     }
  1744. }
  1745.  
  1746.  
  1747. /* MR
  1748.  *--------------------------------------------------------------
  1749.  *
  1750.  * ReconBiMBlock --
  1751.  *
  1752.  *    Reconstructs bidirectionally predicted macroblocks.
  1753.  *
  1754.  * Results:
  1755.  *      None.
  1756.  *
  1757.  * Side effects:
  1758.  *      None.
  1759.  *
  1760.  *--------------------------------------------------------------
  1761.  */
  1762.  
  1763. static void ReconBiMBlock(VidStream *vid_stream, int bnum,
  1764.               int forw_col_start, int forw_row_start,
  1765.               int back_col_start, int back_row_start,
  1766.               int zflag)
  1767. {
  1768.   int row, col, row_size;
  1769.   unsigned char *dest, *past, *future;
  1770.  
  1771.   row_size = vid_stream->mb_width;
  1772.   /* Calculate macroblock row and column from address. */
  1773.   row = vid_stream->mblock.mb_address / row_size;
  1774.   col = vid_stream->mblock.mb_address % row_size;
  1775.  
  1776.   /* If block is luminance block... */
  1777.  
  1778.   if (bnum < 4) {
  1779.     /* Set dest to luminance plane of current pict image. */
  1780.     dest = vid_stream->current->luminance;
  1781.  
  1782.     /* If past frame exists, set past to luminance plane of past frame. */
  1783.     if (vid_stream->past) past = vid_stream->past->luminance;
  1784.  
  1785.     /*If future frame exists, set future to luminance plane of future frame. */
  1786.     if (vid_stream->future) future = vid_stream->future->luminance;
  1787.  
  1788.     /* Establish row size. */
  1789.     row_size <<= 4;
  1790.  
  1791.     /* Calculate row,col of upper left pixel in block. */
  1792.     row <<= 4;
  1793.     col <<= 4;
  1794.     if (bnum > 1) row += 8;
  1795.     if (bnum & 1) col += 8;
  1796.  
  1797.     forw_col_start = (forw_col_start>>1) + col;
  1798.     forw_row_start = (forw_row_start>>1) + row;
  1799.     back_col_start = (back_col_start>>1) + col;
  1800.     back_row_start = (back_row_start>>1) + row;
  1801.  
  1802.   }
  1803.   /* Otherwise, block is NOT luminance block, ... */
  1804.  
  1805.   else {
  1806.     /* Establish row size. */
  1807.     row_size <<= 3;
  1808.  
  1809.     /* Calculate row,col of upper left pixel in block. */
  1810.     row <<= 3;
  1811.     col <<= 3;
  1812.  
  1813.     forw_col_start = (forw_col_start>>2) + col;
  1814.     forw_row_start = (forw_row_start>>2) + row;
  1815.     back_col_start = (back_col_start>>2) + col;
  1816.     back_row_start = (back_row_start>>2) + row;
  1817.  
  1818.     /* If block is Cr block... */
  1819.  
  1820.     if (bnum == 4) {
  1821.       /* Set dest to Cr plane of current pict image. */
  1822.       dest = vid_stream->current->Cr;
  1823.  
  1824.       /* If past frame exists, set past to Cr plane of past image. */
  1825.       if (vid_stream->past) past = vid_stream->past->Cr;
  1826.  
  1827.       /* If future frame exists, set future to Cr plane of future image. */
  1828.       if (vid_stream->future) future = vid_stream->future->Cr;
  1829.     }
  1830.     /* Otherwise, block is Cb block... */
  1831.  
  1832.     else {
  1833.       /* Set dest to Cb plane of current pict image. */
  1834.       dest = vid_stream->current->Cb;
  1835.  
  1836.       /* If past frame exists, set past to Cb plane of past frame. */
  1837.       if (vid_stream->past) past = vid_stream->past->Cb;
  1838.  
  1839.       /* If future frame exists, set future to Cb plane of future frame. */
  1840.       if (vid_stream->future) future = vid_stream->future->Cb;
  1841.     }
  1842.   }
  1843.  
  1844.   dest += (short)row * (short)row_size + col;
  1845.   past += (short)forw_row_start  * (short)row_size + forw_col_start;
  1846.   future += (short)back_row_start * (short)row_size + back_col_start;
  1847.  
  1848.   if (!zflag)
  1849.     BIMcm_reconstruct(dest, past, future, &(vid_stream->block.dct_recon[0][0]), row_size);
  1850.   else
  1851.     BIM_reconstruct(dest, past, future, row_size);
  1852. }
  1853. #endif
  1854.  
  1855.  
  1856. /* MR
  1857.  *--------------------------------------------------------------
  1858.  *
  1859.  * ProcessSkippedPFrameMBlocks --
  1860.  *
  1861.  *    Processes skipped macroblocks in P frames.
  1862.  *
  1863.  * Results:
  1864.  *    Calculates pixel values for luminance, Cr, and Cb planes
  1865.  *      in current pict image for skipped macroblocks.
  1866.  *
  1867.  * Side effects:
  1868.  *    Pixel values in pict image changed.
  1869.  *
  1870.  *--------------------------------------------------------------
  1871.  */
  1872.  
  1873. static void ProcessSkippedPFrameMBlocks(VidStream *vid_stream)
  1874. {
  1875.   int row_size, row, col, addr, off;
  1876.  
  1877.   /* Calculate row sizes for luminance and Cr/Cb macroblock areas. */
  1878.  
  1879.   row_size = vid_stream->mb_width << 4;
  1880.  
  1881.   /* For each skipped macroblock, do... */
  1882.  
  1883.   addr = vid_stream->mblock.past_mb_addr + 1;
  1884.   row = (addr / vid_stream->mb_width)<<4;
  1885.   col = (addr % vid_stream->mb_width)<<4;
  1886.  
  1887.   for (; addr < vid_stream->mblock.mb_address; addr++, col+=16) {
  1888.  
  1889.     /* Calculate macroblock row and col. */
  1890.     if(col>=row_size) col=0, row+=16;
  1891.  
  1892.     /* For each row in macroblock luminance plane... */
  1893.  
  1894.     off = ((short)row * (short)row_size);
  1895.     PMB1_reconstruct(vid_stream->current->luminance + off + col,
  1896.              vid_stream->future->luminance + off + col,
  1897.              row_size >> 2);
  1898.  
  1899.     /* For each row in Cr, and Cb planes... */
  1900.  
  1901.     off >>= 2;
  1902.     off += (col>>1);
  1903.     PMB2_reconstruct(vid_stream->current->Cr + off,
  1904.              vid_stream->current->Cb + off,
  1905.              vid_stream->future->Cr + off,
  1906.              vid_stream->future->Cb + off,
  1907.              row_size >> 3);
  1908.   }
  1909.  
  1910.   vid_stream->mblock.recon_right_for_prev = 0;
  1911.   vid_stream->mblock.recon_down_for_prev = 0;
  1912. }
  1913.  
  1914.  
  1915.  
  1916.  
  1917. /* MR
  1918.  *--------------------------------------------------------------
  1919.  *
  1920.  * ProcessSkippedBFrameMBlocks --
  1921.  *
  1922.  *    Processes skipped macroblocks in B frames.
  1923.  *
  1924.  * Results:
  1925.  *    Calculates pixel values for luminance, Cr, and Cb planes
  1926.  *      in current pict image for skipped macroblocks.
  1927.  *
  1928.  * Side effects:
  1929.  *    Pixel values in pict image changed.
  1930.  *
  1931.  *--------------------------------------------------------------
  1932.  */
  1933.  
  1934. static void ProcessSkippedBFrameMBlocks(VidStream *vid_stream)
  1935. {
  1936.   int row_size, half_row, row, col, ccol, crow;
  1937.   int right_half_for, down_half_for, c_right_half_for, c_down_half_for;
  1938.   int right_half_back, down_half_back, c_right_half_back, c_down_half_back;
  1939.   int addr;
  1940.   int recon_right_for, recon_down_for, recon_right_back, recon_down_back;
  1941.   int c_right_for, c_down_for, c_right_back, c_down_back;
  1942.   unsigned char forw_lum[256];
  1943.   unsigned char forw_cr[64], forw_cb[64];
  1944.   unsigned char back_lum[256], back_cr[64], back_cb[64];
  1945.   int row_incr, half_row_incr;
  1946.   int h;
  1947.   char *lum;
  1948.  
  1949.   /* Calculate row sizes for luminance and Cr/Cb macroblock areas. */
  1950.  
  1951.   row_size = vid_stream->mb_width << 4;
  1952.   half_row = row_size >> 1;
  1953.   row_incr = row_size >> 2;
  1954.   half_row_incr =  half_row >> 2;
  1955.  
  1956.   /* Establish motion vector codes based on full pixel flag. */
  1957.  
  1958.   recon_right_for = vid_stream->mblock.recon_right_for_prev;
  1959.   recon_down_for = vid_stream->mblock.recon_down_for_prev;
  1960.   if (vid_stream->picture.full_pel_forw_vector) {
  1961.     recon_right_for <<= 1;
  1962.     recon_down_for <<= 1;
  1963.   }
  1964.  
  1965.   recon_right_back = vid_stream->mblock.recon_right_back_prev;
  1966.   recon_down_back = vid_stream->mblock.recon_down_back_prev;
  1967.   if (vid_stream->picture.full_pel_back_vector) {
  1968.     recon_right_back <<= 1;
  1969.     recon_down_back <<= 1;
  1970.   }
  1971.  
  1972.  
  1973.   /* If only one motion vector, do display copy, else do full
  1974.      calculation. 
  1975.   */
  1976.   /* Calculate motion vectors. */
  1977.   
  1978.   if (vid_stream->mblock.bpict_past_forw) {
  1979.     right_half_for = recon_right_for & 0x1;
  1980.     down_half_for = recon_down_for & 0x1;
  1981.     recon_right_for >>= 1;
  1982.     recon_down_for >>= 1;
  1983.     c_right_for = recon_right_for >> 1;
  1984.     c_down_for = recon_down_for >> 1;
  1985.     c_right_half_for = recon_right_for & 0x1;
  1986.     c_down_half_for = recon_down_for & 0x1;
  1987.     
  1988.   }
  1989.   if (vid_stream->mblock.bpict_past_back) {
  1990.     right_half_back = recon_right_back & 0x1;
  1991.     down_half_back = recon_down_back & 0x1;    
  1992.     recon_right_back >>= 1;
  1993.     recon_down_back >>= 1;
  1994.     c_right_back = recon_right_back >> 1;
  1995.     c_down_back = recon_down_back >> 1;
  1996.     c_right_half_back = recon_right_back & 0x1;
  1997.     c_down_half_back = recon_down_back & 0x1;
  1998.     
  1999.   }
  2000.   /* For each skipped macroblock, do... */
  2001.   addr = vid_stream->mblock.past_mb_addr + 1;
  2002.   row = (addr / vid_stream->mb_width)<<4;
  2003.   col = (addr % vid_stream->mb_width)<<4;
  2004.  
  2005.   for (; addr < vid_stream->mblock.mb_address; addr++, col+=16) {
  2006.     
  2007.     /* Calculate macroblock row and col. */
  2008.     if(col>=row_size) col=0, row+=16;
  2009.  
  2010.     /* Calculate upper left pixel row,col for luminance plane. */
  2011.     crow = row >> 1;
  2012.     ccol = col >> 1;
  2013.     
  2014.     /* If forward predicted, calculate prediction values. */
  2015.     
  2016.     if (vid_stream->mblock.bpict_past_forw) {
  2017.       
  2018.       ReconSkippedBlock16(vid_stream->past->luminance, forw_lum,
  2019.               row, col, row_size, recon_right_for, recon_down_for,
  2020.               right_half_for, down_half_for);
  2021.       ReconSkippedBlock8(vid_stream->past->Cr, forw_cr, crow, ccol, half_row,
  2022.              c_right_for, c_down_for, c_right_half_for, c_down_half_for);
  2023.       ReconSkippedBlock8(vid_stream->past->Cb, forw_cb, crow, ccol, half_row,
  2024.              c_right_for, c_down_for, c_right_half_for, c_down_half_for);
  2025.     }
  2026.     /* If back predicted, calculate prediction values. */
  2027.     
  2028.     if (vid_stream->mblock.bpict_past_back) {
  2029.       ReconSkippedBlock16(vid_stream->future->luminance, back_lum,
  2030.               row, col, row_size, recon_right_back, recon_down_back,
  2031.               right_half_back, down_half_back);
  2032.       ReconSkippedBlock8(vid_stream->future->Cr, back_cr, crow, ccol, half_row,
  2033.              c_right_back, c_down_back,
  2034.              c_right_half_back, c_down_half_back);
  2035.       ReconSkippedBlock8(vid_stream->future->Cb, back_cb, crow, ccol, half_row,
  2036.              c_right_back, c_down_back,
  2037.              c_right_half_back, c_down_half_back);
  2038.     }
  2039.  
  2040.     h = ((short)crow * (short)half_row) + ccol;
  2041.     lum = (vid_stream->current->luminance + ((short)row * (short)row_size) + col);
  2042.  
  2043.     if (vid_stream->mblock.bpict_past_forw && !vid_stream->mblock.bpict_past_back) {
  2044.  
  2045.       PSB1_reconstruct(lum, forw_lum, row_incr);
  2046.  
  2047.       PSB2_reconstruct(vid_stream->current->Cr + h, vid_stream->current->Cb + h,
  2048.                forw_cr, forw_cb, half_row_incr);
  2049.  
  2050.     } else if (vid_stream->mblock.bpict_past_back && !vid_stream->mblock.bpict_past_forw) {
  2051.  
  2052.       PSB1_reconstruct(lum, back_lum, row_incr);
  2053.  
  2054.       PSB2_reconstruct(vid_stream->current->Cr + h, vid_stream->current->Cb + h,
  2055.                back_cr, back_cb, half_row_incr);
  2056.     } else {
  2057.       PSB3_reconstruct(lum, forw_lum, back_lum, row_size);
  2058.  
  2059.       PSB4_reconstruct(vid_stream->current->Cr + h, vid_stream->current->Cb + h,
  2060.                forw_cr, back_cr, forw_cb, back_cb, half_row);
  2061.     }    
  2062.   }
  2063. }
  2064.  
  2065.  
  2066.  
  2067. /* MR
  2068.  *--------------------------------------------------------------
  2069.  *
  2070.  * ReconSkippedBlock -- width 16
  2071.  *
  2072.  *    Reconstructs predictive block for skipped macroblocks
  2073.  *      in B Frames.
  2074.  *
  2075.  * Results:
  2076.  *    No return values.
  2077.  *
  2078.  * Side effects:
  2079.  *    None.
  2080.  *
  2081.  *--------------------------------------------------------------
  2082.  */
  2083.  
  2084. static
  2085. void ReconSkippedBlock16(unsigned char *source, unsigned char *dest, int row, int col, int row_size,
  2086.              int right, int down, int right_half, int down_half)
  2087. {
  2088.   source += ((short)(row + down) * (short)row_size) + col + right;
  2089.  
  2090.   if ((!right_half) && (!down_half))
  2091.     RSB1_reconstruct(dest, source, row_size);
  2092.   else 
  2093.     RSB2_reconstruct(dest, source, source + right_half + ((short)row_size * (short)down_half), row_size);
  2094. }
  2095.  
  2096. /* MR
  2097.  *--------------------------------------------------------------
  2098.  *
  2099.  * ReconSkippedBlock -- width 8
  2100.  *
  2101.  *    Reconstructs predictive block for skipped macroblocks
  2102.  *      in B Frames.
  2103.  *
  2104.  * Results:
  2105.  *    No return values.
  2106.  *
  2107.  * Side effects:
  2108.  *    None.
  2109.  *
  2110.  *--------------------------------------------------------------
  2111.  */
  2112.  
  2113. static
  2114. void ReconSkippedBlock8(unsigned char *source, unsigned char *dest, int row, int col, int row_size,
  2115.             int right, int down, int right_half, int down_half)
  2116. {
  2117.  
  2118.   source += ((short)(row + down) * (short)row_size) + col + right;
  2119.  
  2120.   if ((!right_half) && (!down_half))
  2121.     RSB3_reconstruct(dest, source, row_size);
  2122.   else 
  2123.     RSB4_reconstruct(dest, source, source + right_half + ((short)row_size * (short)down_half), row_size);
  2124. }
  2125.  
  2126.  
  2127. /*
  2128.  *--------------------------------------------------------------
  2129.  *
  2130.  * DoPictureDisplay --
  2131.  *
  2132.  *    Converts image from Lum, Cr, Cb to colormap space. Puts
  2133.  *      image in lum plane. Updates past and future frame
  2134.  *      pointers. Dithers image. Sends to display mechanism.
  2135.  *
  2136.  * Results:
  2137.  *    Pict image structure locked if displaying or if frame
  2138.  *      is needed as past or future reference.
  2139.  *
  2140.  * Side effects:
  2141.  *    Lum plane pummelled.
  2142.  *
  2143.  *--------------------------------------------------------------
  2144.  */
  2145.  
  2146. static void DoPictureDisplay(VidStream *vid_stream)
  2147. {
  2148.  
  2149.   /* init display, if needed */
  2150.  
  2151.  
  2152.   if(!vid_stream->display_is_initialized)
  2153.   {
  2154.     ResizeDisplay(vid_stream->h_size, vid_stream->v_size);
  2155.     vid_stream->display_is_initialized = TRUE;
  2156.   }
  2157.   /* Convert to colormap space and dither. */
  2158.  
  2159.     if(ditherType == GRAY_DITHER || ditherType == CYBERGFXGRAY_DITHER)
  2160.         vid_stream->current->display = vid_stream->current->luminance;
  2161.     else
  2162.       DoDitherImage(vid_stream->current->luminance, vid_stream->current->Cr,
  2163.                                     vid_stream->current->Cb, vid_stream->current->display,
  2164.                                     vid_stream->mb_height <<4, vid_stream->mb_width <<4);
  2165.  
  2166.   /* Update past and future references if needed. */
  2167.  
  2168.   if ((vid_stream->picture.code_type == I_TYPE) || (vid_stream->picture.code_type == P_TYPE)) {
  2169.     if (vid_stream->future == NULL) {
  2170.       vid_stream->future = vid_stream->current;
  2171.       vid_stream->future->locked |= FUTURE_LOCK;
  2172.     } else {
  2173.       if (vid_stream->past != NULL) {
  2174.     vid_stream->past->locked &= ~PAST_LOCK;
  2175.       }
  2176.       vid_stream->past = vid_stream->future;
  2177.       vid_stream->past->locked &= ~FUTURE_LOCK;
  2178.       vid_stream->past->locked |= PAST_LOCK;
  2179.       vid_stream->future = vid_stream->current;
  2180.       vid_stream->future->locked |= FUTURE_LOCK;
  2181.       vid_stream->current = vid_stream->past;
  2182.       ExecuteDisplay(vid_stream);
  2183.     }
  2184.   } else
  2185.     ExecuteDisplay(vid_stream);
  2186.  
  2187. }
  2188.  
  2189.  
  2190.  
  2191. /*
  2192.  *--------------------------------------------------------------
  2193.  *
  2194.  * ToggleBFlag --
  2195.  *
  2196.  *    Called to set no b frame processing flag.
  2197.  *
  2198.  * Results:
  2199.  *      No_B_Flag flag is toggled from present value to opposite value.
  2200.  *
  2201.  * Side effects:
  2202.  *      None.
  2203.  *
  2204.  *--------------------------------------------------------------
  2205.  */
  2206.  
  2207. void ToggleBFlag(void)
  2208. {
  2209.   No_B_Flag ^= 0xffff;
  2210. }
  2211.  
  2212.  
  2213. /*
  2214.  *--------------------------------------------------------------
  2215.  *
  2216.  * TogglePFlag --
  2217.  *
  2218.  *    Called to set no p frame processing flag.
  2219.  *
  2220.  * Results:
  2221.  *      No_P_Flag flag is toggled from present value to opposite value.
  2222.  *
  2223.  * Side effects:
  2224.  *      None.
  2225.  *
  2226.  *--------------------------------------------------------------
  2227.  */
  2228.  
  2229. void TogglePFlag(void)
  2230. {
  2231.   No_P_Flag ^= 0xffff;
  2232. }
  2233.